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NOISE REGULARIZATION AND COMPUTATIONS FOR THE 
1-DIMENSIONAL STOCHASTIC ALLEN CAHN PROBLEM 

MARKOS A. KATSOULAKIS, GEORGIOS T. KOSSIORIS, AND OMAR LAKKIS 



Abstract. We address the numerical discretization of the Allen-Cahn prob- 
lem with additive white noise in one-dimensional space. Our main focus is to 
understand the behavior of the discretized equation with respect to a small 
"interface thickness" parameter and the noise intensity. The discretization is 
conducted in two stages: (1) regularize the white noise and study the reg- 
ularized problem, (2) approximate the regularized problem. We address (1) 
by introducing a piecewise constant random approximation of the white noise 
with respect to a space-time mesh. We analyze the regularized problem and 
study its relation to both the original problem and the deterministic Allen- 
Cahn problem. Step (2) is then performed leading to a practical Monte-Carlo 
method combined with a Finite Element-Implicit Euler scheme. The resulting 
numerical scheme is tested against theoretical benchmark results concerning 
the behavior of the solution as the interface thickness goes to zero. 

1. Introduction 

Stochastic partial differential equation (SPDE) models arise in numerous ap- 
plications ranging from materials science, surface processes and macromolecular 
dynamics [Coo70, Spo89], to atmosphere and ocean modeling [LN03] and epidemi- 
ology [Dur99]. These models are typically derived from finer and more detailed 
models where unresolved degrees of freedom are represented by suitable stochastic 
forcing terms. There are also some notable rigorous derivations from microscopic 
scales in special asymptotic regimes [BPRS93, MT95, e.g.]. 

An important class of models consists of the stochastic Ginzburg-Landau models 
which are typically obtained from microscopic lattice models for a suitable order pa- 
rameter (e.g., spin), by statistical mechanics renormalization arguments combined 
with detailed balance laws. 

Numerical simulation of these nonlinear SPDE's constitutes an important re- 
search issue. On the practical side, one is interested in having efficient, reliable 
and not too complex numerical codes which can be used either in the context of 
Monte Carlo methods or for sample paths simulations of the physical models of 
phase transition [WB95, KM99, KKOl, ShaOO, e.g.]. From a more theoretical view- 
point, understanding the issues arising from the discretization of SPDE's, in a more 
general setting than phase separation, both through finite difference or finite ele- 
ment schemes, turns out to be a non-obvious departure from numerical schemes for 
deterministic models [ANZ98, Gy699, DZ02, BTZ04, ST03]. 
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In this paper, our focus is on the numerical simulation of the stochastic AUen- 
Cahn problem, which is one of the simplest models exhibiting the phenomena of 
interface formation and nucleation. The stochastic Allen-Cahn problem is an ad- 
hoc white noise perturbation of the deterministic Allen-Cahn, given by 

(1.1) dtu{x,t) - d^^u{x,t) + ,U{u{x,t)) ^ eW^tW{x,t), bTX(ED,te [0, oo) 

where D = (—1, 1) C M^, e > and /e is an odd nonlinearity scaled by and 
dxtW is the space-time white noise (see §2.1 the details). This is a stochastic version 
of the well-known deterministic Allen-Cahn problem describing the evolution in 
time of a polycrystalline material [AC79]. We take boundary conditions of Neumann 
type and the initial condition to be a resolved profile; we refer again to §2 for the 
details. Note that this equation, with white noise, is tractable only in 1 spatial 
dimension, which is the case we will study. In higher space dimension, one has to 
consider noise which is colored in space. 

Equation (1.1) is a type A model in Halperin's classification [HH77]. It is non- 
conservative in the order parameter u and exhibits both nucleation and interface 
formation, whilst retaining a relatively simple structure without multiplicative or 
conservative noise terms encountered in type B models, such as the Cahn-Hilliard- 
Cook equation [KM99] . 



While a thorough discussion of (1.1) is given in §2, it is worth mentioning here 
that this SPDE, with the white noise term, is well-posed only in one space dimen- 
sion. Two important pieces of work concerned with the analytic and probabilistic 
aspects of (1.1) are those of Funaki [Fun95] and Brassesco, De Masi & Presutti 
[BDMP95]. In both papers, the authors study the asymptotic behavior of the solu- 
tion processes as e — 0. In particular, it turns out that, under suitable time-space 
rescaling, the solution with initial value taken to be (roughly speaking) a step func- 
tion, converges (in an appropriate probabilistic sense) to the step function with its 
jump point performing a Brownian motion. 

Though finite difference schemes have been used for simulations [KM99, KKOl], 
we follow here a finite element approach. The reason driving us to understand 
finite element methods (FEM) for such equations is that FEM constitute a quite 
flexible tool, especially for problems in higher dimensions where one may have to 
deal complex geometries. Also, finite elements are naturally suited for adaptive 
schemes where fine scales may be resolved only on small portions of the domain 
in order to obtain a reasonable accuracy. We believe that understanding the FEM 
in a non-adaptive one-dimensional setting will pave the way to more sophisticated 
studies. 

Our strategy to formulate a finite element scheme for (1.1), follows an idea 
introduced for linear problems by Allen, Novosel & Zhang [ANZ98], and consists 
in two steps: 

1. regularize the noise term dxtW, by replacing it with a somewhat smoother 
approximate white noise dxtW; 

2. discretize the regularized problem. 

This approach allows us to conduct a rigorous analysis of the approximation. It 
makes the subsequent finite element discretization straightforward. Note that a 
finite difference variant based on our regularization is also possible. 

Our first task, carried out in §3, is to construct a regularization, denoted dxtW{x, t), 
of dxtW{x,t) (appearing in (1.1)) with respect to an underlying uniform partition, 

X J^p, of the space-time domain D x I. In the spirit of FEM, this regularization 
process consists of a projection of the white noise onto an appropriate space of 
piecewise constant space-time functions, which may be viewed as the mixed deriva- 
tives of hat functions. This idea, which has been successfully used in the context 
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of the linear heat equation [ANZ98], leads to the regularized problem 

(1.2) dtu{x,t) - d^^u{x,t) + f,{u{x,t)) = e^d^tW{x,t), for x e D,t e [0,oo). 

Notice that dxtW is still a stochastic process in space-time, but it is much smoother 
than the white noise which allows equation (1.2) to be interpreted in the usual PDE 
sense pathwise. 

In §4, after recalling some basic properties of problem (1.2) and its solution, 
we prove Theorem 4.4, which states that the solution of the regularized problem 
converges — in an appropriate sense — to the solution of the original SPDE (1.1) as 
the space-time partition becomes infinitely fine. 

Next, in §5, we relate the solution of the regularized problem to the deterministic 
solution of the Allcn-Cahn equation. Our main result here. Theorem 5.4, proved 
for 7 > 3, indicates that the regularization parameters have to be sufficiently small 
for the noise to be captured in the numerical computations. In fact, according to 
this Theorem the weaker the noise, the finer one must take the space-time mesh, 
in order to see the noise effects. This is due to the fact that for a fixed space-time 
mesh and e — >■ 0, the distance between the regularized stochastic solution u and 
the deterministic solution, g, is of higher order in e than the distance between u 
and u. Our proof makes use of the spectrum estimates of the linearized elliptic 
differential operator —dxx + feil), derived independently by Xinfu Chen [Che94] 
and de Mottoni & Schatzman [dMS95] 

We note that while numerical schemes for the stochastic Allen-Cahn involving 
a spectral approach to white noise have been analyzed [Liu03] this is, up to our 
knowledge, a first analysis using projection methods to regularize the white noise. 

Step 2 of our strategy is accomplished in §6, where we derive a simple finite 
element scheme for the regularized problem (1.2). This is a scheme which uses 
piecewise polynomial finite elements to discretize the space variable and an im- 
plicit (backward) Euler scheme to discretize the time variable. Related numerical 
schemes have been thoroughly analyzed and successfully applied in the context of 
the deterministic Allen-Cahn problem [FP03, KNS04, FW05] and for the stochastic 
linear heat diffusion problem [ANZ98]. It is for the first time, up to our knowledge, 
that this scheme is employed in a stochastic and nonlinear setting. The issues of 
regularity of the regularized solution and the convergence of the FEM are objects 
of our current research. 

In §7, we test our scheme in combination with a Monte Carlo simulation. The 
test consists in reconciling the computational results with the theoretical results 
obtained by Funaki [Fun95] and Brassesco, Dc Masi & Presutti [BDMP95] inde- 
pendently. Our benchmarking procedure consists in tracking the so-called center 
of a resolved profile of the Allen-Cahn equation as time evolves, performing sta- 
tistics thereon and comparing them with the probabilistic results coming from the 
theory. The following conclusions are drawn: (1) The robustness of the Monte 
Carlo method depends on the noise intensity, the lower is the noise the higher is 
the observed robustness. (2) The noise has to be resolved satisfactorily in order to 
see stochastic effects. In contrast with the first conclusion, the lower the noise the 
more has one to resolve the mesh in order to see the noise. This is in competition 
with the need to have a fine mesh in order to resolve the transition layer, due to 
the structure of the solution of the Allen-Cahn equation. (3) The behavior cap- 
tured by the numerics is consistent with the theoretical results; in particular, the 
Mueller-Funaki time scale 1 -|- 27 (see 5.6 for the details) and the corresponding 
Brownian motion diffusion coefficient are clearly exhibited by our numerical results. 
We close with some computations that capture the drift of the interface, modeled 
by the Allen-Cahn equation. This drift, typical of the stochastic solution is quite 
fast with respect to the deterministic case where the solutions are metastable states. 
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2. Set up 

2.1. Noisy Allen— Cahn problem. Wc will study an initial-boundary value prob- 
lem associated with the semilinear parabolic partial differential equation with ad- 
ditive white noise, known as the stochastic (or noisy) Allen-Cahn equation given 
by (1.1). The nonlinearity /j is the derivative of an even coercive function with 
exactly two minimum points. A function such as is known as a double-well 
potential and. for sake of conciseness, we focus on the model potential explicitly 
defined by 

(2.1) F,(0 = ^(e'-1)', force M. 

Here e € R"*" is a scaling parameter. The term dxtW is the space-time Gaussian 
white noise, which can be defined as the mixed distributional derivative of a Brow- 
nian sheet W [Wal86, KX95]. The parameter 7 S M models the intensity of the 
white noise and plays a delicate role in the analysis, as e — > 0. 

The presence of the right-hand side makes (1.1) a randomly perturbed version 
of the AUen-Cahn equation which is a stochastic PDF (SPDE). A solution of such 
an equation has to be interpreted in the stochastic sense. That is, for each t, the 
solution u{-,t) is understood as a random process on an underlying probability 
measure space {n,^,P) with values in a suitable function space defined on D. 
Equation (1.1), supplemented with the initial condition 

(2.2) u{x,0) = uo{x),VxeD, 
and with the Neumann boundary conditions 

(2.3) a^u(-l,t) = a^w(l,t) = VtGM+, 

defines the stochastic Allen-Cahn problem. For simplicity, we assume that the 
initial condition uq is smooth enough and satisfies the boundary conditions. In §5 
we shall focus on a more particular class of initial conditions known as resolved 
profiles. 

2.2. Space-time stochastic integral. One can give a mathematically rigorous 
definition of a solution of the stochastic Allen-Cahn problem (l.l),(2.2)-(2.3) as a 
distribution- valued process [Wal86, KX95]. However, we find it more convenient, 
as in the case of the white noise generated from a Brownian motion, to work with 
the stochastic integral with respect to the Brownian sheet W denoted by "J • dW'^ 
[Wal86, §11] [KX95, Ch. 3]. In our doing so, we bear in mind the formal relationship 

(2.4) / / f{x,t)dxtW{x,t)dxdt^ / / f{x,t)dW{x,t) 

Jo J D Jo J D 

that will inspire the weak formulation (2.10) and the definitions in §3. In the 
particular case where / is the characteristic function of a Borel-measurable set 
A e ^(M+ X D) of Lebesgue measure |A| < 00 the following basic property of the 
stochastic integral is satisfied 



(2.5) / / dW{x,t)^W{A)e^{Q,\A\), 

Jo J A 

i.e., M^(A) is a Gaussian random variable with mean zero and variance \A\} 

Since we are interested in numerical solutions, we consider the time domain to 
be a bounded interval / = [0,r], for some fixed T € R+. A fundamental property 



"'^For /J G R, o" 6 M+ we denote by N (/x, cr'^) the class of normally distributed (or Gaussian) 
random variables of mean and variance on the space fl. 
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of the stochastic integral is the foUowing weU-known h2-isometry, which holds for 
the Ito integral, 

2" 



(2.6) 



f{x,t)dW{x,t) 



I JD 



= E 



for any ^/^-mcasurablc / G L2(/ x Z? x fi), where 

(2.7) = a {W{A) ■ Ae.^ilx D)} , 

is the sigma-field (or sigma-algebra) generated by W up to time t, and E denotes 
the expectation with respect to (O, P).^ 

A useful consequence of (2.6) is that 
(2.8) 



/ JD 



f{x,t)dW{x,t) / / g{y,s)dW{y,s) 



I JD 



f{x, t)g(x, t) dx dt 



I JD 



for any -measurable /, g e L2(/ x Z) x il). In the special case where / and g are, 
respectively, the characteristic functions of two Borel sets A and B E x D), 
with \A\ , \B\ < oo, (2.8) implies 

(2.9) Coy{W{A),W{B)) ^\AnB\. 

2.3. Integral solutions. By multiplying (1.1) with a test function g C^{D x 
(0, oo)) and using the formal relation (2.4), one can write the problem in the usual 
weak form'^ 



(2.10) 



{udt<l> - dxudx(f> - fe{u)(j>) + 



^dW = 0. 



Despite the above formulation being quite useful, especially for studying a numerical 
scheme, it is not very convenient to nail down the concept of solution. A rather 
more convenient way to give rigorous meaning to (1.1) is to look for an integral 
solution of an equivalent integral equation [DPZ92, Doe87, FJL82, Wal86], as we 
briefly illustrate next. 

Introduce first the corresponding boundary value problem for the stochastic lin- 
ear heat equation [DPZ92, Wal86] 

dtZ - dx,xZ = d^tW, in £» X M+ 

(2.11) Z{x,0) = 0, onD 

9,Z(l,t)-9,Z(-l,t) = 0, Vie [0,oo). 

The solution to this problem can be defined as the Gaussian process in space-time 
produced by the stochastic integral 

ft 



(2.12) 



Zt{x)=Z{x,t) 



Gt^,{x,y)dW{y,s), 



JD 



where G is the heat kernel for the corresponding homogeneous Neumann problem. 
In our one-dimensional particular case, G can be explicitly written as 

. , 7rfc(a- + l) nkiy + l) -ir^kH 
(2.13) Gt{x, y) = 4^(2 - 5„^) cos cos ^y-^ exp 

fe=0 

where Sq is the Kronecker symbol. 



In compliance with the standard practice in stochastic differential equations, we write explic- 
itly the probability variable oj g f2 as an argument to random variables only when necessary in 
order to avoid confusion. 

■^Whenever the meaning is clear from the context, for sake of conciseness, we often drop the 
variables "x,t" and, in non-stochastic integrals, also the corresponding elementary terms "d". 
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The integral solution of (1.1) can then be defined as a solution of the equivalent 
integral equation 

(2.14) 

u{x,t) ^ - / Gt-s{x,y)ff,{u{y,s))dyds+ Gt{x,y)uoiy) dy + e"^ Zt{x). 
Jo JD JD 

It is known that such a solution exists uniquely as a C*'(£')-valued continuous 
process, t i— > u{.,t), adapted to Zt, provided the initial condition uq satisfies the 
Neumann boundary conditions [BDMP95, Wal86, FJL82]. In this article we use this 
concept of solution which we refer to simply as the solution of Problem (1.1),(2.2)- 
(2.3) and we will denote it by u. Notice that u is also referred to by some authors 
as the Ginzhurg- Landau process [BDMP95]. 

For the aims set in this paper, namely, in order to study the error of convergence 
of an approximation of the solution of (1.1), we will need a uniform bound for u. 
While in the deterministic case such a bound is direct consequence of the maximum 
principle, in the stochastic case one cannot expect to have a uniform bound in the 
whole probability space. However, a bound on a set with large probability controlled 
by e will suffice for our needs. We present an extension of a previously known result 
of Brassesco et al. [BDMP95, Pro.5.2]. 

2.4. Lemma (Probabilistic maximum principle) Let 7 > —1/2. For each T > 
and Kq > there exist ci, C2, (5o > such that if ||uo||l < 1 + <5o then 

(2.15) p\ sup \\uit)\\^ > l + A'ol <ciexp(-C2/ei+''^). 

[te[o,T] J 

Proof We reduce the proof to that of [BDMP95, Pro.5.2] by introducing the time- 
space rescaling: t ^ t/e^ and x ^ x/ (v^e) and extending the solution periodically 
to the whole space as to obtain the proper barrier function. Since we are dealing 
with the more general case 7 > —1/2, while they deal with the case 7 = only, 
we retrace the salient points of their proof. The barrier function v satisfies the 
following equation — corresponding to [BDMP95, (5.12)]: 

(2.16) 5..4a... + 2..-3.^-.3^2-VV-/2a.,W^. 
Consider now the function 

(2.17) V{x,t) = f\xp{-2{t-s))Htf\x,y)dWiy,s), 

Jo 

where Hj:^^'' is the Green operator defined by 

(2.18) exp(-2t)ijtf ^ = (^dt - ^a,. + 2Id 

with homogeneous boundary conditions on (— l/(A/2e), l/{\/2e)). By using equa- 
tion [BDMP95, (5.2)] with A = exp(— (7 + 1/2)) and adapting properly the proof 
of [BDMP95, Lemma 2.1] we can easily conclude that for each 6 > there exist ci 
and C2 > such that 

(2.19) PU sup e'^+^/'^V{x,t) > b\ <cicxp{-c2/e^+'^-^). 

1 t<Te-2,a;GK I 



The rest of the proof is now standard. 
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3. White noise approximation 

In order to introduce a finite element method (FEM) that approximates a solu- 
tion of (2.14), we first need to to obtain a weak formulation in the standard sense 
of PDE and FEM. This is not possible with the presence of the white noise, so we 
regularize first the problem by replacing the white noise with a smoother stochastic 
term. Our technique is inspired by that of Allen, Novosel k, Zhang [ANZ98] for the 
linear heat equation. 

3.1. A piecewise constant approximation of the white noise. Consider a 
tensor-product partition of the space-time domain, S!^ x ^p, where cr, p (z K+ and 

{Dm ■ Dm (a;,„_i , x^), TO G [1 : M]} , 

and J^p := {/„ : /„ := [t„_i,t„), n £ [1 : N]} , 

are, respectively, a space-domain, and a time-domain, partition; each one of these 
partitions is uniform, that is 

(3.2) Xjn — x„i-i ~ a, Vto. G [1 : M] and t„ — t„_i = p, Vn S [1 : N] 

and xo = ~1, xm = 1, = and In = T. We denote by Xm = 1_d„ and 
ifn = l/„ the characteristic functions of the space subdomains and time subdomains 
respectively. 

The (piecewise constant) approximation of white noise, abbreviated by AWN 
below, is given by the random space-time function 

N N 

(3.3) d,tWix,t) = Y,Yl 

{x)(pn (t) 

n—1 m—1 

where the coefficients are the random variables defined by 

(3.4) j7m,„:=— / / Xm{x)^n{t)dW{x,t). 

o-p J I Jd 

In the sequel we will use the shorthand 

(3.5) / [ f{x,s)AW{x,s)^ I I f{x,s)d,tW{x,s)dxds, 
Jo Jd Jo Jd 

in spite of the integral being taken in the classical, non-stochastic, sense. 

3.2. Lemma (Moments and independence of the AWN coefficients) TJie coefRcients 
fjm.n defined in (3.4) are i.i.d. N(0, l/{ap)) variables. 

Proof From the definitions of firn,n and property (2.5) we have 

^m,n = — I I Xni{x)(pn{t)dW{x,t) 

J I Jd 

(3.6) ^^11 dW(M) = 2^^i^ 

<^P Jin Jd„, 

V p ) \ op 

To show independence compute the covariances for to, to' £ [1 : M\ and n, n' G 
[1 : TV], using (2.8), as follows 

(3.7) 



X7nfndW / Xm-'Pu'dW 

iiJd JiJd 



Xm{x)Xm' ix)ipn{t)ipn' (t) dx dt 
O^'On'^^P, 



1 1 JD 

cm cn 
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where Sj is the Kronecker symboL I 
The AWN satisfies two important technical properties that we state and prove 
next. 

3.3. Lemma (Approximate Ito-type inequahty) For all deterministic functions f € 
L2(/ X D) the following holds true 



(3.8) 



E 



fix,t)dWix,t) 



I J D 



< J J f{x,tydxdt. 



Proof Lemma 3.2 and some manipulations yield 



E 



fix,t)dWix,t) 



I J D 



E 



JlJD ^„ 



T\m,nXra(x)(l)n{t) dx dt 



f{x, t) dx dt 

Yl ( / / / ) + 2 Yl ^™^"^m'« 



/„ -ID 



iL -Id 



< 



E 



per \Jj^ Jd 



f{x,tfdx dt. 



I JD 



In the next-to-last step we use the Cauchy-Schwarz inequality. 
3.4. Remark Lemma 3.3 and (2.6) imply that 



(3.9) 



E 



f{x,t)dW{x,t) 



I JD 



< E 



fix,t)dW{x,t) 



I JD 



In other words, the L2-type regularity properties of the AWN will be, at the worse, 
the same as those of the white noise itself. 

Since we will need bounds on space-time norms of the AWN, but in probabilities 
rather than in expectation, we establish the following basic result. 

3.5. Lemma (Loo(L2) and L2(L2) bounds for the AWN) For each A' > we have 
(3.10) 



P\ sup \\d,tW{t)\\ <K\> 



1-- ( 1 
P 



2 \ l/<^-l 



exp 



K' 



-P 



and 



(3.11) 

Proof We proceed in steps. 



~2 



exp 



2 
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Step 1. Recall that AI = 2/a and N ^ T / p. By the definition of dxtW we have, 
for each t e [0, T] and n e[l : N] such that t e In, that 



(3.12) 



\d..tW{t)\ 



1.2(D) 



M ^ M 

'yy-t 1 ' itr-in 1 



Vm.n ' 



m— 1 m— 1 

where the rj^^n G N (0, 1). In order to conclude, we will obtain a condition on the 
right-hand side that makes it smaller than K^, for all n € [1 : iV]. 
Step 2. For each n g [1 : iV] we consider the random variable 

AI 

(3.13) F„ E ^" 



Im , n • 



Notice that, in view of Lemma 3.2 for n n' , Hn and Hn' are independent. Let us 
fix n for a while and find an event for which Hn < pK^. By Lemma 3.2 and a basic 
probability fact [Bil95, Pbm. 20.16], the random variable Hn has a chi-squarcd 
distribution with M degrees of freedom. Its density is given by 

z*^/2-iexp(-z/2) 



(3.14) 



for z > 0, 



2M/2r(M/2) 

and for z < 0, where F is the Euler Gamma- function. Thus we have 

rPK' 



(3.15) 



P{Hn<pK^} 



1 



2M/2r(Af/2) 7o 



z^^/2-icxp(-z/2)dz. 



Step 3. Wc prove next a lower bound on this integral in the case where M is even, 
the odd case being similar. Let y play the role of pK^ and consider for each k G Nq 
the integral 

rv 

exp(-z/2) dz. 



(3.16) 

An integration by parts yields the recursive expression 
(3.17) h = Ikh-i - 2/ exp(-y/2), 

which allows, by an inductive argument, to see that 

fc! 



(3.18) 



/fc = 2'^'+U-!-2E 



(fc-z)l 



/-'2^exp(-y/2). 



An easy manipulation with the binomial formula implies that 

(3.19) h > 2''+^kl (1 - (1 + j//2)'= exp(-y/2)) . 

Taking fc — M/2 — 1 in the above and recalling the definition of Ik and (3.14) it 
follows that 

(3.20) p{Hn< pK^}>l- + exp 
which implies 



(3.21) 



P{Hn<pK^} > 



I- 1 



pK' 



\ M/2-1 



exp 



pK^ 



Step 4. To conclude the proof, we introduce the event 

N 

(3.22) - fl {Hn < pK^} , 
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and we observe that, in view of (3.12), on ilj^ we have 



(3.23) 



Vie [0,r]. 



On the other hand, using the independence of iJ„, n e [1 : A^], the simple fact that 
(1— ^)^>1 — for ^ < 1 and (3.21) we can estimate the probabihty 

N 

Pi^l)='[[P{Hn<pK^} 



(3.24) 



> 



> 



1-1 



.2X M/2-1 



exp 



1 -iV 1 + 



2\ M/2-1 



exp 



pK^ 
pK^ 



N 



By replacing N T/p and M = 2/cr we get (3.10). 

Step 5. Estimate (3.11) is obtained simply by using (3.21) with pK^ and M re- 
placed by K'^ and MN respectively. 

I 

3.6. Remark (alternative proof) As pointed out by one of the referees, it is possible 
to prove Lemma 3.5 more directly, by using martingale inequalities. 

3.7. Remark (interpretation of (3.10) and (3.11)) We may rewrite the term ap- 
pearing in (3.10), as 



(3.25) 



T f 



l/n-l 



exp 



TcxpF{p, a,K). 



A practical way to use such a result is by fixing first T, p,a € R"*" and then requiring 
a big enough K such that T exp F{p,a, K) <C 0. This is made possible by the 
fact that ImiK^oo F{p, a, K) = — oo for any fixed p,(T & R^. The same type of 
observation is valid also for the (3.11). 

4. The regularized solution 

We now introduce the regularized solution to problem (1.1), (2.2)-(2.3), which 
we obtain by replacing the white noise by the AWN in (1.1). The role of the 
regularized problem is pivotal in devising a numerical scheme to approximate the 
stochastic Allen-Cahn problem. We discuss the approximation properties of this 
regularization with respect to the original problem. 

4.1. Definition (regularized solution) The regularized solution^ u, of the noisy 
Allen-Cahn problem is the unique continuous solution of the integral equation 



(4.1) u{x,t) 



Gt-six,y)fe{u{y,s))dy ds 



JD 



Gt{x,y)ua{y)Ay + e^ / Gt.,{x,y)dW{y,s). 

D JO JD 

4.2. Lemma (maximum principle for regularized solutions) For fixed T, Kq > 0, 
there exist So, ci, C2 > 0, independent of e, sucli tliat if |1uo||l (_d) !i 1 + <5o tlien 



(4.2) 



P { sup \\u{t)\\^ < 1 + Ao ^ > 1 - ci exp(-C2/ei+27). 
te[o,T] 
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Proof We follow exactly the proof of Lemma 2.4, by observing that (3.9) ensmes 
that all the estimates for the stochastic integrals of the white noise can be "trans- 
lated" in corresponding estimates for the integrals of the approximate white noise. 
The constants appearing in this Theorem can be therefore taken to be the same 
that appear in §2.4. I 

4.3. Remark (regularized solution is strong solution) Notice that the regularized 
solution u of (4.1) is in fact a weak solution in the PDE sense, i.e., u{t;Lj) G 
H^(Z?) and dtu{t;uj) e L2(-D) for all t S (0,r] and ui € H,, and the following weak 
formulation is satisfied: 

{dtu{t;uj),(f>) + {d^u{t;uj),d^(l)) + {f,{u{t;uj)),(f>) 

(4.3) =e'^{d.,tWit;Lu),^), V0 e HJ(7^), i G (0, T], 
and u{0; uj) = uq, 

for each lu £ il (the notation (•, •) indicating the inner product in Li2{D)). Indeed, 
each one of the AWN's realizations, dxtW{u!), is a picccwise constant space-time 
function. For each such realization the usual regularity theory for scmilincar para- 
bolic equations with picccwise continuous data can be applied and the correspond- 
ing weak formulation written down [LSU68] . 

Our next goal is to show that the regularized approximate solution converges to 
the solution u. For this we will estimate the regularization error 

(4.4) e{x,t) = u{x,t) - u{x,t), 

in terms of the white noise regularization parameters a and p, and show that it 
converges to zero in an appropriate sense. 

4.4. Theorem (convergence to the stochastic solution) For a fixed T, tiiere exist 
constants ci, C2, Ci and C2 such that for each e € (0, 1) there correspond an event 

and a constant > such that 

(4.5) P{n^) > 1 - 2ci exp(-C2/e^+2'') and 



(4.6) / ( / f \u-u\A dP<cJc,p'^^ + C2^],ya,p>0 



e{x,t)^l I Gt-s{x,y){fe{u{y,s))~ f^{u{y,s)))dyds 
(4.7) ^ 

+6-^1 f Gt^.,ix,y){dWiy,s)~ dWiy.s)) 



pi/2 

Proof We proceed by steps. 

Step 1. By the integral representations of u, (2.14), and u, (4.1), we can represent 
the error as an integral too: 



ID 

for all {x,t) £ D X (0,T]. So our task now is to bound the terms in the left-hand 
side of (4.7) in the appropriate norm. 

Step 2. In view of the maximum principle for both the exact solution, §2.4, and 
the approximate solution, §4.2, there exists an event ilf C such that 

(4.8) P{nf) > 1 - 2ci exp(-C2/ei+'^) 
and 

(4.9) C {Wum^^^D) - < 3, Vt g [0,T]} . 
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The choice of the number 3 is quite arbitrary here. In fact any number of the form 
1 + Kq with Kq > will do, with the appropriate change of constants. This and 
the local Lipschitz continuity of / imply that 

28 

(4.10) \Mu) - /,(u)| <-\u-u\,o^ nf. 

Step 3. Working now on the event and introducing the functions 



(4.11) £(r) ~ / / e{x,tydxdt 

Jq Jd 

(4.12) (/.(r) := 



JD 



JD 



Gt-s{x,y){dW{y,s)- dW{y,s)) dxdt 

I 

for all r e [0,T], wc infer from (4.7) that 

(4.13) eir)<2 J^(^J^^\Gt-,ix,y)\'^e{y,s)dyds^ dx dt + 2e^'' ^{r) . 

The integral in (4.13) can be bounded, using the Cauchy-Schwarz inequality, by 
(4.14) 

2^/ f (f [ \Gt^s{x,y)\' dyds [ [ e{y , dy ds) dx dt ^ [ z{t)e{t)dt 
£ Jo JD \JQ JD Jq Jd / Jo 

where 

(4.15) z{t):^2^[ f f \Gt^s{x,y)\^ dydsdx. 

£ JdJo Jd 

Inequality (4.13) implies 

(4.16) e(r) < 0(r) + /" z{t)e{t)dt, 

JO 

for each r ^ I . Applying the Gronwall lemma to this inequality we obtain 

(4.17) e(T)<exp( / z{t) dt] e"^ <t^{T) < C,<t>{T), 



e(T) < exp ^ 



where — by estimating the heat kernel — the constant is given by 



(4.18) C, := e^-* 



e ' exp 



V 12 £4 



Step 4. By summing with respect to P on the event both members of this 
inequality we obtain 

(4.19) [ [ [ \u-uf dxdtdP <C, [ (j){T)dP <G,E[(j){T)]. 
Jfi"^ Jo J d Jn^ 

We conclude by observing [ANZ98, Lem. 2.3] that there exist Ci, C2 > 0, depending 
only on T, such that 

2 

(4.20) E[cj){T)] < Cip^/^ + C2^. 
Thus we established that 

(4.21) X i / 'i^'i^'i-P^'^^ (^i'°'^' + ^2^ 
as we claimed. 
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4.5. Remark (About the constant C^) Theorem 4.4 insures that, for fixed T and e, 
the approximate solution u converges to u as p, c — > 0. The constant Cc appearing 
in the estimate depends exponentially on both 1/e* and T, thus for small e, or large 
T, this might force us to take very small p and a. This fact should to be taken into 
account in practice. The bound we have proved seems to be pessimistic though, as 
the choice of a and p, used in our subsequent numerical experiments, indicates. 

4.6. Remark (convergence rate) Recalling that p and a can be thought as dis- 
cretization (in addition to regularization) parameters, the convergence rate found in 
(4.6) is in acccordance with standard results for linear parabolic equations [ANZ98, 
e.g.]. Note that the need for p = Ctr^, the so-called "parabolic space-time scaling", 
we obtain the right balance between the two terms appearing on the right-hand 
side of (4.6). 



In this section we focus on the relation between the regularized stochastic AUen- 
Cahn problem (4.3) and the deterministic version. The reason to do this, is to find, 
in an analytical setting, what conditions should be taken on the regularization 
parameters, p and u for the noise to be captured in the regularized equation. 

We mainly show that the the error between u and the deterministic solution to 
the Allen-Cahn problem, g, in an appropriate probability-Loo (0, T; L2) sense, is of 
order O(e^) as e for fixed 7 > 3 and p,a > 0. In §5.6 we give an interpretation 
of this result as an evaluation of the risk of obtaining a poor resolution of the 
noise for fixed p and a and too small e. This poor resolution may lead to the 
disappearance of the stochastic effects in the regularized equation, even after the 
appropriate rescaling, because u becomes much closer to q than u, with respect to 
e. This point is further investigated numerically in §7. 

Our proof makes use of a spectrum estimate result for the linearized Allen-Cahn 
operator [Che94, dMS95], which is recalled in Theorem 5.2 and the L2{D) estimate 
on the noise given by 3.5. The proof's technique is a continuous data dependence 
result for parabolic equations based on a Bernoulli-Gronwall type argument, in the 
spirit of Feng & Wu [FW05]. 

The result holds for 7 > 3 and it is an open problem, as far as we know, to find 
the critical 7 for which the result ceases to hold. 

5.1. Deterministic solution and resolved profiles. Denote by q be the (clas- 
sical) solution of the problem 



We consider also the function of space only go defined as the unique solution to 



(In fact it is qo = tanh.) We will assume from now on that uq is a resolved profile 
solution^ which is defined to be an e-lincar perturbation of an e-rescalcd and shifted 
go- That is, for all x S uq{x) = qo{{x — xo)/e) + epo{x) where xq € Z?, and 
Po is such that uq satisfies the Neumann boundary conditions Uq(±1) ~ 0. With 
this choice of initial condition the linearization of the operator u —dxxU + /e(u) 
about q enjoys the following spectral property. 



5. The regularized solution's limit for e — > 



(5.1) 
(5.2) 
(5.3) 



dtq - dxxq + fe{q) = 0, in D x / 
q(0) = Uq, on D 
5,g(t,0) = 9:,g(t,l) = 0, tel. 



(5.4) 



% + fM = in M, go (±00) = ±1 and go(0) = 0. 
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5.2. Theorem (Spectrum estimate [Che94, dMS95]) There exists a constant Aq > 
such that for any e E (0, 1] we have 

(5.5) Wd.nliD) + ifM)^,^) > -Ao Ml(D) : V0 e R\D). 

It is also a well-known consequence of the maximum-principle that if |uo| < 1 
(which is the case when uq is a resolved profile) then \q\ < 1. 
The main result of this section is 

5.3. Lemma (continuous dependence for the regularized-deterministic error) There 
exists a bounded and non- increasing function Ki : [0, oo) — ^ M and a constant K2, 
both depending only on Aq, such that 

(5.6) 11^2(0 - g(i)!lL,p) <i^2e' 
provided 



(5.7) / ||a,tVK(s)||' exp(-(3 + 2Xo)s) ds < Kiit)e''-^\ 



lL2(r>) 
forte [0,T]. 

Proof We divide the proof in several steps and we denote in it ||-||l2(-D) simply by 

IMI- 

Step 1. We start by deriving an energy inequality for the error 

(5.8) e :^ u — q. 

Since u satisfies the weak formulation (4.3) and g is a classical solution, we can 
write the following PDE in its weak formulation for e: 

(5.9) {dte, 0) + {d^e, 8,4,) + {f'M% 4>) - {d.tW , 0) - ^ (fe^ 0) , V</) e Yi\D), 
where 

(5.10) f ■.= iq + e^2q + u. 
Testing with e in (5.9) we obtain 

(5.11) (a*e,e) + ||a,e||V(/;(g)e,e) < {d^^W ,e) ~ ^ {f,e^) . 

Step 2. The next step is to bound the terms in the right-hand side of (5.11). The 
first term can be written as 

(5.12) e^d,tW{t),e{t)) <^^\\d.tW{t)\\^ + ■ 

To produce a bound on the second term of the right-hand side of (5.11) we use 
(5.10), valid in 1 spatial dimension, to obtain 

(5.13) {f{t),e{tf) = 3 {q, e{tf) + MmliD) ■ 

By the fact that \q\ < 1 and the Sobolev embedding ll^{D) ^ LodD), valid for 
Z) C M, the first term on the right hand side can be bounded using 

|3(g,e-(i)')| < 3||e-(t)||L^(^) \\e{t)f < Ci WmikHD) Wmf 

(5.14) /-2 

<>^ie'\\m\\lHD) + ^\\m\f 

where Ci is 3 times the Sobolev embedding constant for D and Ai := min{l,Ao} 
(the reason for this choice will be apparent in the next step) with Ao from (5.5). 
As a consequence we have 



(5.15) — 2(^'^^)^Aie 



1 / ^\ . , 2 11 -|i2 , C2 II _||4 1 II -||4 
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where C2 = CiV(4Ai). 

Step 3. Owing to the spectrum estimate (5.5) and the fact that f'{q) > —1 we 
have 

Wd.ef + {f'M), e') =:A = il- e')A + e^A 

>-{l-e')Xo\\ef + e'\\d,e\f~\\er 

(5.16) = - ((1 - e')Ao + 1) + \\d,ef 

= -{l + Xo)\\ef+e^{\\d,ef + Xo\\ef) 
> -(l + Ao)||e||VAie2j|g||^,^^^. 

The inequahties (5.11), (5.12), (5.16) and (5.15) lead to 

(5.17) id,||e-(t)||^-(l + Ao)||e(t)f 



< V P^*^{t)f + I \\e{t)f + ^ ||e(t)r - 1 WemliD) Vt £ /• 



Consider, for the rest of the proof, the following notation: 
(5.18) 

g{t) := \\e{t)\\\ a (3+2Ao), 6 2C2/e^ r(t) 6^7 ||5^^vK(t)f -||e(i)||^^(^) 
Then (5.17) implies 

(5.19) g'{t) < ag{t) + bg{tf + r, Vi G [0, T]. 

Step 4. To proceed we will apply a Bernoulli differential inequality technique, 
which generalizes the Gronwall Lemma, in order to get a bound on g{t). We follow 
Feng & Wu [FW05, Lem.2.1] 

Fix a i € [0,T] and let 



(5.20) q{s) := / exp(— ar)r(r) dr and p(s) -.^ pt{s) = (g(t) — gi(s)) exp(as), 

Jo 

for all s G [0,t]. Since 

(5.21) p'{s) = d,,pt(s) = -r(s) + ap{s) and p{s) > 0, for s e [0,t], 
we may write 

(5.22) d,(g(s)+p(s))<a(g(s)+p(s))+5(.g(s)+p(s))2, forse[0,t]. 
Introducing z{s) := l/{g{s) +p{s)), we can rewrite this inequality as 

(5.23) z'(s) +az(s) > -6, Vs e [0,i]. 
Multiplying by exp(as) and integrating over [0,t\ we obtain 

(5.24) zit) > z{0) exp(-ai) - ^1 - exp(-«t)) ^ 

a 

Noting that .g(0) = |le(0)|| = 0, pt{0) = g{t) and pt(t) = 0, this yields 

J_ > a - &g(^)(cxp(at) - 1) 
g(t) ~ aexp{at)g{t) 
We now invert both sides of this inequality, under the sufficient condition that 

(5.26) a - bg{t){cxp{at) - 1) > 0, 
and we get 

(5.27) \m\\'< «'^-p(«^)^^w 



a — bg(t)(cxTp{at) — 1) 
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Step 5. To conclude we want to interpret more explicitly this result. Let us replace 
first (5.26) by the sufficient condition 

(5.28) a-bg{t){exp{at)-l)>d{t), 

for some S{t) > that will be chosen appropriately. This is equivalent to 

(5.29) eM-as)ris) ds (= git)) < ^(Jp^jj^l^) - 
This can be ensured if we assume 

(5.30) / ||9,tiy(s)f exp(-as)ds < ^""^^^^^^^^ 



(5.33) S{t) 



II - V /II ' - 2C2(exp(at) - 1) 

Under this condition we obtain the bound 

<^-) ii^'""'^ 2£,k") -' 

Step 6. We conclude by taking 

(5.32) J(t) := max |a — exp(ai) + 1, ^1 , 

i.e., 

a — exp{at) + 1, for i < 

for t>ta 

where ta = log(l + a/2) /a. Then, after putting 

(5.34) K,{t) n"Ml,«/2(exp(aO-l)} ^ ^ 

condition (5.29) may be replaced by 

(5.35) / ||a:rtVK(s)f cxp(-as)ds < A"i(t)e^-2T, 
Jo 

a condition under which we have, from (5.31) 

(5.36) ||e(t)||' <i^2'e^ 
where := supjQ q^-j Ki{t)^ exp(at) < 00. 



As a consequence of this estimate we state the following result, which, roughly 
speaking, implies that in order for the noise to have the chance of persisting in the 
limit, as e — >■ 0, the parameters p, a must also go to zero. 

5.4. Theorem (low-intensity approximate white noise) There exists a constant C = 
C(Ao) such that for all fixed j > 3, p,a,T > we have 

(5.37) liui P \ Lo e n : sup ||?l(-; uj) - q\\^^ < Ce^ I = 1. 



[0,T] 



Proof Let C = /i2 in (5.3). Choosing to use an L2(0, t; L2(£'))-norm estimate 
for the white noise, for the estimate (5.6), in view of the monotonicity of Ki, it is 
enough to assume the sufficient condition 

nt 

(5.38) / \\da:tW{s)\\^ ds < KiiT)e^-^''. 
Jo 

According to (3.11), this condition is satisfied with probability 

(5.39) l-(l + ^e^~A'^'''''\xJ-'^e^-^'' 
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Since 7 > 3, 6 — 27 < and, for fixed cr, p and T > 0, it is possible to make this 
probability arbitrarily close to 1 for e > small enough, as claimed. I 

5.5. Remark Note that, it is also possible to obtain a variant of 5.4 by employ- 
ing (3.10) instead of (3.11). This leads to slightly better lower estimates of the 
probability for the same e > for longer time T. 

5.6. White noise resolution by the AWN. We now describe one interpretation 
Theorem 5.4. 

Note firstly that, in the theorem's statement, the AWN regularization parameters 
a and p are kept fixed while e — > 0. It is well known that, if uq is a resolved profile 
with center at (see §7.2 for a definition of center), i.e., mq = tanh((x — a;o)/2e), 
as e — > 0, the solution q = q'^ oi (5.1) converges, in an appropriate sense, to the 
stationary step function Xxq ■= l(i;o,oo) ~ l(-oo,2;o)- Theorem 5.4 is saying that 
for fixed p, cr > and for e — > 0, the solution, u, of the regularized problem (4.3) 
converges to this stationary step function. 

Second we note that, owing to a result by Funaki [Fun95, Theorem 8.1] or a 
similar one by Brassesco et al. [BDMP95], there exists a stochastic process {t,uj) i-> 
^ti^) such that 

(5.40) limPiwefi: sup ||u(-, t; w) - xc.^M > -5 ^ 0, 

for each fixed 6 > 0, where Xxo is the step function defined above, and (i, w) i— >■ is 
a stochastic process which converges as e — > 0, in an appropriate sense (in law), to 
the standard Brownian motion rescaled as to have diffusion coefficient yJcQf^l'^^'^ 
where cq = 3^2/4. Of course, as e — > 0, (w) — >■ a;o, where xq is the center of the 
initial condition mq; this implies that the limit of u and u are consistent, as e — > 0, 
even when p and a are kept fixed. 

Suppose now that one wishes, in view of Theorem 4.4, to use the regularized 
solution u(t\ instead of u(t), to approximate the diffusion coefficient of the process 
t I— > One way of doing this would be to approximate (numerically) u(t,aj), for 
w € f2 (or a discrete analog), find its center, if it exists, and finally compute 

its average (excluding solutions that have no center) and its variance over G fi. 
The resulting variance, rescaled appropriately, i.e., t i— > Var[^(]/(e^"*'^'''), for the 
computation to be meaningful, one should see, asymptotically as e -> 0, a linear 
function t H> c^t. This rescaling, which we shall call the Mueller-Funaki rescaling 
[Fun95] , is necessary in order to get a result that is essentially independent of e and 
thus easy to visualize. 

Theorem 5.4 tells us that, for a fixed p, cr > 0, the rate of convergence of u — > g is 

0(e'^/2). Since the distance IE mit) — Xi\ ) j as can be seen using a piecewise 

constant approximation of tanh, is 0(e^/^), it follows that u is closer to q than 
Xii and that any statistics conducted on u may lead to wrong results. This is a 
strong indication, which is confirmed by the numerical results in §7, that in order to 
capture the stochastic effects the parameters u and p must be chosen as functions 
of e. Note that a similar conclusion can be derived from Theorem 4.4 in case the 
dependence of proves to be effective, but the nature of this similar conclusion 
has its roots in deterministic considerations rather than stochastic ones. 

Although we have proved Theorem 5.4 for values of 7 > 3, it is natural to expect 
similar results for lower values of 7. In fact, our numerical experiments in 7 indicate 
that this is the case. 
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6. An Euler-Galerkin finite element scheme 

We introduce now the finite clement discretization of tlie regularized problem 
(4.3). 

6.1. Discretization partitions. We begin by introducing the space and time par- 
titions 

■■= {D'„, : D'„^ := {x'„^.„x'J, m G [1 : M']} , 

and A: := {K, ■ I'n ■= n e [1 : N']} . ' 

These partitions do not necessarily coincide with the partitions and J^p used 
for the regularization procedure in §3.1. Bearing in mind that this setting could 
be further generalized, we limit ourselves here to the case where the numerical dis- 
cretization partitions, and J^k, are refinements of the white noise regularization 
partitions and J^p, respectively. For each D'^ G S'h there exists Di G i^o- such 
that D'^ C Di etc; this determines a unique mapping /i : [0 : M'] [0 : M], such 
that D'j^ C -Dp(m)- For simplicity, we also assume that the partitions are uniform 
and that the meshsize and timestep are denoted respectively by h and k. The rea- 
son we do not make these partitions coincide is that for the finite element method's 
convergence analysis it may prove useful to have more involved couplings of the 
type h = h{a) and k = k{p). In this article we consider only the simplest situation 
possible where h = a and k = p. 

6.2. Finite element space and the discrete scheme. Let V C HJ(D) be the 
space of continuous piccewise linear functions associated with the partition S>ii, 
we define the (spatial) semi-discrete solution as the time-dependent random finite 
element function L/ : [0, T] x f2 V which solves the SDE 

(6.2) {dtU{t), V) + {d,u{t),d.,v)+{uu{t)), V) - v),yv^Y,te [o, r]. 

We discretize further this SDE in the time variable by taking a semi-implicit Euler 

scheme in time associated to the partition / = {io} U Um ^"i 

(6.3) 

,i^U(9.f/",a.F)+(/,(t/"-i),F) = (a,til',F),VFG V, tG [o,T]. 

The scheme is called semi-implicit, in that the linear part is treated implicitly while 
it is explicit in the nonlinearity. This means that at each timestep only one linear 
problem has to be solved and no nonlinear solver is needed. 

In practice, it is more practical to use a modified version of (6.3) given by 

^g^^ ^ 1 ,v^ + {d.u-,d,v) + {nxu-')u\v) 

= (/.'(t/"- !)[/"-! - /,([/"-!), y) + {d.,tW, v^) , vy G V, t G [0, t]. 

which allows bigger timesteps k [KNS04] . Note that this amounts to a linearization 
involving one step of the Newton method to solve the nonlinear (fully implicit) 
backward Euler scheme. 

6.3. The linear time-stepping system. Let us indicate the basis functions of 
V by for m G [0 : A/']; that is the piecewise linear continuous function such 
that $„^(a;^) = JJ", for Z G [0 : M']. If we indicate by = {yC) the vector of 
nodal values corresponding to the discrete solution [/" at time i", that is U^{x) = 
'Y^^=o'^Tn^rn{x), then, we can translate (6.4) in the following matrix form 



(6.5) 



m=0 ' 



-M + A + ^N(u''- 
k €^ 



k 
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where M, A are the usual finite element mass and stiffness matrices, respectively, 
7V(m"^^) and ^(m"^^) are a "nonlinear" mass matrix and load vector, respectively, 
and w = (wm) a random load vector generated at each time-step. A short calcula- 
tion shows that for m an internal degree of freedom (node) we have 

(6-6) = ^-^(77^(,„)_l + 77M(m)) 

where ^ is the mapping introduced earlier in this section and 77; is a N (0, 1) random 
number for I £ [0 : M], or zero for I = — 1,M -|- 1 (the boundary cases). If the 
partitions S>a and coincide, which will be the case in the next section, then 
h ~ a and (6.6) simplifies to 

(6.7) Wrn = + ^™)- 

This is the form that we employ in our computations below. From now on, we will 
consider .^p and J^fc to coincide, i.e., k = p in all our computations. 

7. Computations 

We discuss now numerical simulations of the Allen-Cahn problem, using the 
scheme (6.5). Our main purpose is to match the behavior of the exact stochastic 
solution with the behavior of the Monte-Carlo type numerical solution. By "exact 
behavior" we refer to the theoretical results of Funaki [Fun95] and Brassesco et 
al [BDMP95]. Our numerical experiments are directed towards testing relations 
between the various parameters e, 7 and p, a. In addition to taking k = p, we also 
take h = (7. 

7.1. Monte Carlo simulations. We use the finite element scheme (6.4) in combi- 
nation with Monte Carlo type simulations. For each choice of parameters e, 7 and 
meshsize h, we choose the timestep k ^ and compute between 1000 and 2500 
samples paths, each with a different seed for the pseudo-random number generator. 

Each sample path runs from time to a final time T. At the beginning of each 
run, the random number generator is seeded and the subsequent rjm appearing in 
(6.7) are chosen according to this seed for all the run. The seed for each sample path 
is determined by the clock of the machine at the start of each run (these are also 
recorded for rerunning purposes) . In this section we denote the numerical solution 
(which tacitly depends on e, 7, h, etc.) by (C^")ne[0:7V]i where n corresponds to the 
timestep and w is a discrete sample, i.e., the choice of the initial seed. We indicate 
by fl the discrete sample space, which can be thought of being the set of all initial 
seeds employed. 

7.2. Phase-separation interfaces and centers. Our benchmarking procedure 
consists in comparing the behavior of the center of the discrete solution with that 
of the exact solution. A function v € }1^{D) is said to have a center if v{x) = is 
uniquely solvable with a sign change for x £ D = ( — 1, 1); the solution of v(x) = 
is called center of v. For example the function uo{x) ~ tanh((a; — 2;o)/'\/2e) has 
center xq. 

The center represents a 0-dimensional interface separating two phases in the 1- 
dimensional Allen-Cahn model here considered. Due to the noise a solution with 
center can nucleate, i.e., give rise to new zeroes, during the evolution; in this case 
it makes no longer sense to speak of "the" center and we simply say interface 
position. Such solutions, as we shall see, must be treated carefully in the statistics. 
We note also that an interface may disappear, either by exiting the domain D or 
by annihilating another interface. 
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(a) Interface average (top), E[S"], and its 
rescaled variance (bottom), Var[S"]/20e^"'"^^, 
as functions of time t„, for e = 0.04. Only one 
line, corresponding to the coarsest refinement 
level I = 6, i.e., meshsize h = 2~^, is visible 
as for all the other levels no admissible paths 
survive up to time T = 20 with a center (i.e., 
only one zero). 
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(c) Interface average (top), E[S"], and its 
rescaled variance (bottom), Var[S"]/20e-'^"'"^^, 
as functions of time t„, for e = 0.01. In each 
plot here, since e is even smaller, we see four 
graphs corresponding to all four refinement lev- 
els tested, I = 6, . . . , 9. Note that the vari- 
ance for i = 6 is, practically zero, which is far 
from the exact results. We believe, that this is 
due to two concurrent factors: (1) the meshsize 
is too coarse to resolve the interface layer, (2) 
the meshsize is too coarse to resolve properly 
the noise and the numerical solution is closer 
to the exact deterministic solution rather than 
the stochastic one. Note also that the variance 
time-dependence, which is initially linear as ex- 
pected, starts degenerating at about T = 4. 
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(b) Interface average (top), E[S"], and its 
rescaled variance (bottom), Var[S"]/20e'^+^'', 
as functions of time tn, for t = 0.02. Compar- 
ing this with the case e = 0.004 (Figure 1(a)) 
we see that more graphs arc visible in each di- 
agram, corresponding to the refinement levels 
I = 6, 7, 8. This is due the fact that for smaller 
value of e more sample paths the noise affects 
the solution less. 



-^0 




(d) These are the results from a different per- 
spective, with a two-fold purpose: (1) to check 
the scaling power 1 -|- 27 in the diffusion coef- 
ficient, and (2) to see how the quality of re- 
sults depends behaves with respect to time. 
Each of these graphs show logVar[S"] ver- 
sus loge for the finest refinement level ( = 9. 
Each graph corresponds to a fixed choice of 
the discrete time index, n, among a sequence 
ni < n2 < ■ ■ ■ < rt/, with rt; = 10ni_i. A 
reference line of slope 1 + 27 which is the slope 
that is expected for the other lines, the code to 
pick up. We see that this is reflected by the 
computations, but with the quality worsening 
as time grows. 



Figure 1. Numerical results for 7 = 0.0 



7.3. A benchmark. According to known results [Fun95, Theorem 8.1][BDMP95] 
and the discussion done in 5.6, we expect the center of fJ" € H^(D), if any, to 
perform a (discrete) Brownian Motion, modulo perturbations of order 0(e) and the 
numerical discretization error. Therefore, we will declare our numerical scheme to 
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be acceptable, if we see evidence of this Brownian motion in our computational 
results. 

Based on this observation, our benchmarking procedure consists in tracking the 
center (SJ^)o<n<T/fe C (—1,1) of the computed sample path {U^)o<n<T/ki where 
w ranges in a discrete sample space, and we perform statistics on (5")o<„<T/fe by 
averaging over w. In order to obtain meaningful statistics, we keep the number 
of interfaces constrained to 1. Paths that maintain a center up to the final time 
are called admissible sample paths. If new interfaces are created, or the center 
exits from (—1, 1), during the computation of one sample path, it is rejected as not 
admissible, the computation stopped, and the computation of a new sample path 
is started. The average and variance of the interface position are computed over 
the admissible sample paths. The resulting average interface position, E[S"], and 
its variance, Var[S"] = E[(S")^] — (E[S"])^, are real valued functions of (discrete) 
time tn- 

According to (5.40) and the subsequent remarks, the average position E[S"] must 
be close to and its variance must be close to coe^'^'^'^tn, and this is what we will 
be after in the next section. 

Note that our use of the benchmark is no so strict than the one usually used in the 
context numerical schemes for deterministic nonlinear equations. In the latter case, 
an exact solution is usually readily obtained and the benchmark procedure simply 
consists in measuring the error between the two solutions. In the linear stochastic 
case, a similar approach can be used with the moments of the SPDE's solution 
[ANZ98]. We should stress, however, that in our case, which is both stochastic and 
nonlinear, the analytic knowledge of the exact solution is too scarce, making such 
simple benchmarking nearly impossible. 

7.4. Simulations and results. We run a scries of Monte Carlo tests with various 
combination of parameters as following: 

S = Q 

while 5* < S'inax do 
for Lj =clock time do 

seed of the random number generator with uj 



for 7 = 0.0, 0.2, 0.5 do 

for e = 0.08, 0.04, 0.02, 0.01 do 

for ^ = 6, 7, 8, 9 {I is the refinement level} do 
let h = 2i-' and k ^ 
for n e [1 : T/k] do 
solve (6.5) for 

find EZ such that U^i^Z) = 
if S" exists and is unique then 



else 

break and skip to the next I 
end if 
end for 

declare sample path successful: 



let E[S"]„e« 
let E[(S")2] 



new " 



new 



= mE-r]S + EZ)/S 



S = S + 1 

E[S"] = E[5"]„ew 
E[(S")2] = E[(S")2] 



now 



end for 
end for 
end for 
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end for 
end while 

The results are reported in Figures 1, 2 and 3 for the values of 7 = 0.0, 0.2, 0.5, 
respectively. In each figure the sub- figures (a), (b) and (c), each of which is split 
into a top and bottom part, show the graph, as (discrete) functions of tn G [0,T], 
of the (discrete) average position E[S"] in the top part and its variance VarS" := 
E[(e;")2] - E[E;"]2 rescalcd by l/ei+^T in the bottom part for the values of e = 
0.04, 0.02, 0.01, respectively. (For easier visualization we plot the piecewise linear 
interpolation of discrete functions.) Different lines in each of these 3 diagrams 
correspond to different values of the mesh refinement level I . The absence of a line 
means that the total number of successful sample paths is and no statistics are 
produced. The diagrams (a), (b), (c) correspond, on each figure, to the values of 
e = 0.04, 0.02, 0.01. The plot (d) is designed so to see how well the scheme captures 
the CQC^^'^^t behavior of the variance; each line corresponds to a chosen but fixed 
value of t, and is a log-log plot of the variance for the finest refinement level (Z = 9) 
against e; different lines correspond to different times between and T = 20. For 
comparison we plot the line with slope 1 -I- 27, which represents the scaling power; 
lines parallel to this line mean that the code picks up the right scaling. 

We summarize next the observations we have drawn from the computational 
results that were described. 

Interface motion. The motion of the numerical interface, for sufficient transition 
layer and noise resolution, has the properties predicted by the analytical results. 
Indeed, starting from a resolved profile centered at zero, as seen in Figures 1-3 the 
average position is near zero, whereas the variance, which is expected to be a linear 
function of time, behaves in accordance to the expectations, at least for some initial 
times; the 7, e dependence of the diffusion coefficient cqe^'^'^'^ is clearly captured by 
the numerics, as seen by the diagram (d) of each of the figures. 

Noise resolution. It is well known that for simulating the deterministic Allen-Cahn 
equation with any type of mesh/grid-based schemes, the meshsize has to be smaller 
than e in order to resolve satisfactorily the transition layer about the interface. 
Intuitively, this is due to the fact that the transition layer being of width 0(e), and 
the numerical discretization parameters must be smaller than this width as to have 
a proper resolution of this transition layer. 

The effect of this is seen in each of Figures 1 to 3, as e decreases, the level 
of refinement has to be taken bigger and bigger in order to obtain meaningful 
calculations. 

In the stochastic case, the situation is complicated even more by the noise. 
Indeed, according to Theorem 5.4 and the discussion in §5.6, the discretization 
parameters, in this case h = a and k = p, must be taken small enough as to ensure 
that the noise effects are not lost for small e. Roughly speaking, the discretization 
parameters must be small, not only to resolve the interface, but also to resolve 
the noise and pick up the diffusion of the Brownian Motion. This dependence, 
which is indicated analytically for 7 > 3 by Theorem 5.4, is also reflected in our 
computations for 7 < 3. Indeed, a comparison between Figure 1 and Figure 2 
shows how the level of refinement / = 7 leads to meaningful results for 7 = 0.0 
and all values of e (at least for short times), whereas the same refinement level, for 
the same values of e, but with 7 = 0.2 is insufficient. For 7 = 0.5 in Figure 3 this 
phenomenon becomes yet more apparent. 

We note that computations with 7 > 3 (not shown here) require an extremely 
fine mesh, and thus a very small timestep in view of Remark 4.6, in order to capture 
any of the noise effects. Otherwise the deterministic solution will be computed. In 
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this case, even choices of h that resolve the interface satisfactorily are not enough 
to resolve the noise. The lesson we learn from this, is that the interplay between 
the noise and the nonlinearity can be quite delicate, and not obviously predicted 
from deterministic considerations, in problems such as the stochastic Allen-Cahn 
equation. 

Deterioration of simulations with t big and e small. As observed in the previous 
paragraph, the computed variance depends linearly on time, as expected, but only 
for some initial time. The smaller e, the shorter this time is. This is seen in 
the bottom part of the sub-figures (a), (b) and (c) by the graph's earlier or later 
departure from an initial linear behavior. The computations deteriorate faster for 
small e, e.g., e = 0.0, than they do for bigger e, e.g., e = 0.5. The example with 
e = 0.2 shows an intermediate behavior. 

7.5. Interface drift. To conclude, we add some results of computations, for short 
times, with initial value a resolved profile centered away from 0. In this case, the 
SDE describing the motion of the interface has also a drift term, which drives the 
interface towards the closest boundary of the domain. This drift is clearly seen for 
various choices of the parameters in the top diagrams of Figure 4, where we plot the 
center's position average against time. Short times must be taken, for the statistics 
to make sense, otherwise solutions with centers that exit (or nucleate) cease to 
counterbalance those who stay in. Indeed, in the diagram (c), where the noise 
intensity is quite strong, after an initial drift towards the boundary, the average 
inverts its route and moves away from the boundary. This is due to the fact that 
the statistics become too biased; the SE diagram shows the samples survival with 
respect to time (that is, a discrete-probability space version of the exit-time inverse 
function). The samples survival for the top diagrams is 100% for the statistics in 
the top diagrams. It is worth mentioning that similar observations, using stochastic 
ODE's were made by Shardlow [ShaOO]. 
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(a) Interface average (top), E[S"], and its 
rescaled variance (bottom), Var[S"']/20e^"''^^, 
as functions of time t„, for e = 0.04. Compared 
to Figure 1(a) we see here that, due to the more 
modest intensity of the noise, three refinement 
levels produce enough admissible sample paths. 




(c) Interface average (top), E[H"], and its 
rescaled variance (bottom), Var[S"]/20e^"'"^^, 
as functions of time t„, for e = 0.01. We see 
that the linear dependence on time is clearly 
visible in the variance, but for low refinement 
levels (coarse meshsize) the slope may not be 
the proper one. Also for very low refinement, 
the stochastic dynamics arc not picked up at 
all. 
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(b) Interface average (top), E[S"], and its 
rescaled variance (bottom), Var[S"]/20e^+'^''', 
as functions of time t„, for e = 0.02. All four 
refinement levels are represented here. A clear 
linear trend in the computed variance, up to 
time 10 (with a minor deterioration between 
10 and 20) is visible. This linear dependence is 
the expected behavior from the theory saying 
that Var[H"] coe^+^^tn. 
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(d) Each of these graphs show log Var[H"] ver- 
sus loge for the finest refinement level ( = 9. 
Each graph corresponds to a fixed choice of 
the discrete time index, n, among a sequence 
ni < n2 < . . . < n/, with rt; = 10ni_i. A 
reference line of slope 1 + 27 which is the slope 
that is expected for the other lines, the code to 
pick up. Comparing with the case e = 0.0, an 
improvement in the quality of the behavior as 
tn grows is visible, but still a clear deteriora- 
tion in time is present. 



Figure 2. Numerical results for 7 = 0.2. A comparison with 
Figure 1 is due and leads to two important observations. (1) The 
statistics are more robust in time (especially for e small) in that 
the variance shows a clear linear dependence on t in Figures 2(b) 
and 2(c) and the correct scaling behavior in Figure 2(d). This is 
due to the lower intensity of the noise because 7 is bigger. (2) On 
the other hand, we must work harder with the refinement level in 
order to pick up the stochastic effects. We believe that this is a 
practical aspect of the discussion in §5.6. This poor performance 
for bigger meshsizes could not directly related to the interface layer 
resolution (a deterministic effect) as for the same values of e and 
h but lower values of 7 we obtain meaningful statistics, as shown 
in Figure 1. 
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(a) Interface average (top), E[S"], and its 
rescaled variance (bottom), Var[H"]/20€"'^^"^^, 
as functions of time t„, for e = 0.04. 
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(b) Interface average (top), E[S"], and its 
rescaled variance (bottom), Var[H"]/20e^^''^'', 
as functions of time t„, for e = 0.02. 




(c) Interface average (top), E[H"], and its 
rescaled variance (bottom), Var[S"']/20€^^"^^, 
as functions of time t„, for e = 0.01. The 
linear time dependence of the (rescaled) vari- 
ance is clear here. Note that for level of mesh 
refinement / = 6, 7 the computed solution is ba- 
sically deterministic (average is and variance 
is 0). This is a further indication of the prac- 
tical importance of noise resolution, following 
the observations in fj5.6. 




(d) Each of these graphs show log Var[H"] ver- 
sus loge for the finest refinement level ( = 9. 
Each graph corresponds to a fixed choice of 
the discrete time index, n, among a sequence 
n\ < <,...< ni, with = 10ni_i. A 
reference line of slope 1 + 27 which is the slope 
that is expected for the other lines, the code to 
pick up. The value of 7 is quite high now as to 
keep the statistics robust with respect to time. 



Figure 3. Numerical results for 7 = 0.5. The same comments 
made in Figure 2 apply here, but with an even clearer linear de- 
pendence of the variance upon time and a clearer need for high 
mesh refinement in order to resolve the white noise. Also the com- 
putations arc more robust with respect to lower values of 7 and 
they show almost no deterioration in time. 
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(a) The average of the interface position (center) for 7 = 0.0. An upward drift, that 
is a drift of the interface towards the closest boundary point, is visible in this series of 
computations. The number of the Monte Carlo samples is 1623. Only 3 of these paths 
violate the uniqueness of center condition and are therefore excluded from the statistics 
at the time when that happens. 
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(b) The average of the interface position for 7 = 0.2. Compared with Figure 4(a), an 
outward drift is also detectable here but, due to the weaker noise, intensity it is slower. 
The number of Monte Carlo samples here is 1752, all of which make it up to time T = 2. 
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(c) On the left, wc plot the average of the interface position for 7 = —0.1. In this case, most 
sample paths do not make it with the condition of one interface, or center, up to time T = 2. 
This fact is visualized on the right, where the curve represents the number of "admissible" 
sample paths (i.e., those that have a center) versus time; as time increases the number of the 
admissible sample paths drops from an initial 2500 to almost 1500. This high drop in the 
number of samples makes the statistics unreliable, and explains why the graph on the left 
exhibits no clear outward drift as a result. 



Figure 4. Monte Carlo computations of the drift effect for the 
initial interface position shifted away from zero. All these compu- 
tations are performed for e = 0.01, meshsize h = and timestep 
k = h? . See the text and Shardlow [ShaOO], for example, for details 
about the drift. 



